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Abstract 

We study a 2-dimensional model of fermi fields ip that is closely related to the Gross- 
Neveu model, and show that to leading order in jj a (tptp) condensate forms. This effect 
is independent of the chemical potential, a peculiarity that we expect to be specific to 2 
dimensions. We also expect the condensate to be unstable against corrections at higher 
orders in We compute the Green's functions associated with the composite i/ji/j, and 
show that the fermion acquires a Majorana mass proportional to the gap, and that a 
massless Goldstone pole appears. 



Recently several papers have appeared [1-3] dealing with the properties of QCD at high 
density. The basic procedure is to approximate QCD by a direct four-quark interaction 
term, justifying this either by appeal to instanton effects or to one-gluon exchange. In the 
presence of a chemical potential, this theory admits a condensate of quark-quark pairs, very 
similar to the Cooper pairs that are well-known in the BCS theory of superconductivity. 
This phenomenon, which has been dubbed "color superconductivity," may or may not be 
accessible to experiments on heavy-ion collisions that will be performed over the next few 
years. 

In this paper we shall examine similar phenomena in the context of a one-plus-one 
dimensional model that is a close relative to the Gross-Neveu (GN) model [4]. Prom its 
inception, it has been recognized that the GN model exhibits many of the same features as 
QCD, such as asymptotic freedom and spontaneous chiral symmetry breaking. Moreover, 
the four-fermi interaction is the model - there is no need to regard it as an approximation to 
an underlying gauge theory. Unlike in higher dimensions, in two dimensions this interaction 
is renormalizable, and we shall find, just as in the usual GN model, that coupling constant 
renormalization removes all the divergences that we shall encounter. Another advantage 
is that by judiciously introducing a flavor index i, 1 < i < N, one can insure that the 
mean-field approximation that is commonly used in analyzing the condensate in QCD is 
in the case of the 2-dimensional model justified as the leading contribution in powers of 
1/N. 

There are, however, a couple of peculiarities associated with two dimensions that make 
this GN-like model qualitatively different from QCD. The first is the Coleman-Mermin- 
Wagner theorem [5], which forbids spontaneous symmetry breaking of a continuous sym- 
metry. Whereas in the original GN model the broken symmetry is a discrete one, and 
hence not in conflict with the theorem, in this case the formation of a (ipip) condensate 
breaks fermion number, a continuous symmetry. The same problem arises in the chiral GN 
model [4], where the symmetry is continuous, and in a variety of other two-dimensional 
models where the spontaneous breaking of a continuous symmetry is predicted in leading 
order in 1/N. This means that instabilities must arise in higher order that vitiate the 
prediction of a condensate. However, as Witten has pointed out [6], the 1/N expansion 
may still be an excellent guide to the physics of the model, except for the formation of 
the condensate (what happens in these models is that the condensate "almost" forms, in 
the sense that the pair-pair correlation function decays in the infrared only like a power 
instead of exponentially, and the power vanishes as N — > oo). 



The second peculiarity is, as we shall show below, the chemical potential has nothing 
to do with the formation of the condensate. In higher dimensions the chemical potential is 
crucial, because it gives rise to the Fermi surface at which the gap equation has an infrared 
singularity as the gap goes to zero. It is this feature that insures that the gap equation 
will have a solution for arbitrarily weak coupling. In two dimensions, however, the Fermi 
surface has dimension zero, and the infrared singularity exists whether or not there is a 
chemical potential. In fact, the gap equation turns out to be completely independent of 
the chemical potential. This behavior will be exhibited explicitly below. 

The model we consider is defined by the following Lagrangian: 

£ = ^(Oj ^(0 + 2^ 2 -0« 75 V' (i) V' (l) 75V' (j) • (1) 

ipW is a two-component spinor with a flavor index that takes on N values. Repeated 
flavor indices in eqn. (1) are summed. Because of the unconventional arrangement of 
flavor indices in the second term, the model does not have SU(N) symmetry, but it does 
possess O(N) symmetry, ip^ — > d^ip^' where t? lJ is a real N x N orthogonal matrix. 
The Lagrangian (1) also has a (7(1) symmetry, which we shall find is broken by a (ipip) 
condensate, whereas the O(N) symmetry is kept intact. 

Our representation for the 7-matrices is: 7 = o"i; 7 1 = —icri\ 75 = cr 3 , and it is then 
easy to check that 

^W W ¥ 4) 7^ W) = -l(e af 3^% {i %^M j) ) ■ (2) 

Following the usual Hubbard- Stratonovich procedure, in the form introduced by Coleman 
[7], we add to £ the term 

- ^(Bt - g 2 e a ^i (<) V#°)(5 + 9 2 ^HP) (3) 

which does not affect the physics because B and are simply auxiliary fields. We then 
have 

£ = ^(0 (i p _ ^/)^) _ ±-B^B + B e a/ ^t(^t« _ B t ea ^)^f (4) 



where we have also introduced a chemical potential fx. In anticipation of taking the large 
iV limit, N — > oo with A = g 2 N fixed, we rewrite the second term as — ^B^B. 

Thus the classical (or tree-level) term in V e ff(B^B) is of order N. As we perform a 
perturbation expansion, we observe that additional factors of N arise in 2 ways: (i) the 
B — B^ propagator is proportional to 1/N; and (ii) each closed fermion loop gives a factor 
of N from summing on the flavor index. If we examine the computation of higher-loop 
corrections to the effective potential (i.e. the summation over all one-particle-irreducible 
diagrams with zero-momentum B and B^ external legs) we see that the one-loop term is of 
order N (one fermion loop and no B — B^ propagator) whereas anything else is of order N° 
or lower. Hence the leading contribution is just what we get from keeping the tree and one- 
loop graphs. Moreover, in a path-integral approach, one first integrates out the fermions. 
Then, because the exponent is proportional to N, one can employ the stationary phase 
approximation in the integral over B and B^ to evaluate the integrand at the solution of 
the equations 

dB = dW=°- (5) 

The task of integrating out the fermions is complicated slightly by the ipip and i/j^ip^ terms. 
We observe that 



/ 



Dl j, e MM4>) = det^M , (6) 



where M, is any anti-symmetric matrix. We write the fermion part of our Lagrangian as 

£ = I (^AtP - ij>A T ifi) + V t HV t + VSty (7) 

where 



A = (id + ia 3 d x - /u) a p5 lJ ; (8a) 

A T = (-ido - ia 3 d x - iu) aP 6^ ; (86) 

B = Be a ^ = iB(a 2 ) a(3 S ij ; (8c) 

fit = -iB\a 2 ) a p8 ij . (8d) 



Now we perform a translation, 



ip = x + (9«) 
= x + ^ot T (96) 

where a = \{B^)~ X A T ; this factorizes the x an d path integrals into the product of 
two path integrals, and after some manipulations and discarding an overall factor that is 
independent of B and B\ we obtain 

e O' Bt) = detail + 4A- 1 B(A T )- 1 B^] (10) 

where L^ff denotes the one-loop contribution to the effective action. 
Note that the matrix A -1 is given by 

A - 1( , f d 2 k [fcp-^fcx + ^^e^^) 

[X,V) J (2tt) 2 (ko + fi + iesgnko) 2 -kl 1 J 

where k ■ x = kox° + k±x , and the z e prescription is introduced in the proper way to 
take account of the chemical potential. (A T )~ 1 is the same expression but with x and y 
interchanged. 

To obtain the effective potential, we take B and B' to be constant. It is then con- 
venient to write everything in momentum space, and after some algebra, using r e // = 
—V e ff f d 2 x, we obtain 

NB^B iN f d 2 k 



NB^B iN f 

v '"-— + tJ 



trln 



2iW t , d2) 

, AB^B 

[(k + fi + i e sgnk ) + fcicr 3 ] [(k - /i + i e sgnk ) - kia 3 ] 



where the trace is over the spinor indices only. Setting B^B = M/4 and k = 4A, and 
observing that V e ff(M = 0) = 0, we can write 

1 1 r M dV 

V eff = 4- dM'Z-, (13) 



N Veff = NJo aW1 dW> 



where 



1 dV _ 1 i 
N dM ~ k + 2 



(2^ 



M — /cq + (fei<7 3 +ri 2 -»e 



(14) 



Here the trace is just summation over 0-3 = ±1. 

We note that this integral is logarithmically divergent. We shall deal with this by 
renormalizing k; but first we shall do the ko integral. Let M + (k±as + fx) 2 = uj 2 with 
oj > 0. We have 



/: 



dko 



1 



ko — uj 2 + ie 



dko 



1 



(ko — uj + ie) (ko + uj — ie) 



so 



711 
UJ 



(15) 



1 dV _ 1 

NdM ~ k 



1 

8tt 



dk\tr(- 



1 



(16) 



We renormalize by requiring that at 11 = 0, 



1 d 2 K// , 4 
AT dBdB^ \ B ' B = M °/ 4 - KR 



(17) 



which is the same as 



1 



(18) 



The solution to this is 



l_ = 1 _ 1 r- ^ + sx 

K K R 4n J_ 00 UJo 



(19) 



where uJq = k\ + M , and where for our choice of renormalization prescription, SX 
Some other choice of prescription would yield a different pure number for SX. 



j_ 

4-7T ' 



We then have 



1 dV 1 



1 



N dM kr 8tt 



pOO 


"1 1 " 


/ dk\tr 




I — oo 





+ 5X 



(20) 



The gap equation is just the statement that ^ vanishes: 



Kr 8n 



/•OO 


"1 


1 " 


/ dk\tr 






I — oo 




UJ _ 



(21) 



Note that tr - 



+ 



w v /M+(fc 1 + At ) 2 v /M+(fc 1 - At ) 2 



which is even in fci. Therefore 



Kr 4lT 



where J is the integral 



J 



-f 

Jo 



dk 



+ 



y/M + (k + ^) 2 ^/M + (k-^) 2 VM + k 2 



(22) 
(23) 



If we evaluate J, we can find the particular M that obeys eqn. (22), thereby solving the 
gap equation. We can also obtain the more general expression 



N dM kr 4tt 



as a function of M, and integrate it to obtain V(M) via eqn. (13). 
We find, after some mild computational exertions, that 



J(M) = -In M/M 



and 



(24) 



-jyV(M) = (J- + 5X)M + ^-(InM/Mo - 1) + 0(1) 



(25) 



As advertised, these expressions are independent of /U. The solution to the gap equation 
for our choice of 5X is 



M = M e (1 - 47T/KR) 



(26) 



There is no critical lower bound to kr below which no solution exists. We see from eqn. 
(23) that this is due to the fact that as M — > 0, the expression for J diverges logarithmically. 
This infrared singularity is present in 2 dimensions independent of the value of fi. In higher 
dimensions, we expect this singularity to be present at the Fermi surface, k =| [i |, and to 
disappear as \i — > 0. 

We see from eqn. (26) that as M is increased for fixed M, kr becomes smaller. This 
is an indication that the coupling kr is asymptotically free, just as in the original GN 
model. In fact, it is not hard to show from the renormalization condition (19) with cutoff 
A that the beta function has the form 



-K 2 



2tt 



so that the gap, eq. (26), obeys (2M ^ - (3(k r )^-)M = 0. 

As a consequence of the Coleman-Mermin- Wagner theorem we expect that terms 
which are higher order in jj will destabilize the leading order result, i.e. will give rise to 
contributions that dominate the ones we have found for sufficiently large M . 

We may also [4], [8] compute the Green's functions associated with the fields B and 
. To do this, we need the effective action, not just the effective potential. This can be 
read off from equation (10): 

r e// = J d 4 x(^^)B\x)B(x) - l -Trln[l + 4A~ 1 BA~ 1 B J( ] . (27) 

Here we have defined A = o 2 A 7 02- Because the gap equation is independent of fx, we 
shall simplify our task somewhat by setting /j, = 0; then A T = —A, and 

A = id + ia 3 d! (28) 
A = -id + ia 3 d! . (29) 

Note that A A = AA = d$ - d\ . 

To obtain the desired Green's functions, we perform the following sequence of steps: 

(a) We write 



B = B + B' 



(30a) 



=B + B't 



(306) 



where, as above, So is a solution of the gap equation (M = AB 2 ). 

(b) We expand r e // to second order in B' and B'^ . The linear terms in B' and S't will 
cancel because of the gap equation. The coefficients of the quadratic terms will be the 
inverses of the Green's functions that we seek. 

(c) We observe that by introducing the real and imaginary parts of B' : B' = <f>i + ifa, 
B'^ = (pi — i(f>2, the off-diagonal terms will disappear; i.e. there will be no mixed 0102 
terms. In fact, what we find is: 

r e/ /^-4iv J <Px<P y {Mx)Mv) J -0j2 e ~ ip<x ~ y) ^-^+^ 

+Mx)Mv) J 0j2 e ~ ip ' {x ~ y) ll - ^-(p)]> ( 31 ) 



where 

« (No/ MPO + Kg) - Kijpi + «i) ± M] 

±[P) ~ J (2nri-Kl + Ki + M][-(K +Po) 2 + (K 1 +p i r + M] 



(32) 



Now 5>±(p) is logarithmically divergent, but so is ^, and using the gap equation it is easy 
to see that the divergence cancels, along with all residual dependence on kr and 5X. 

(d) The integrals defining &±(j>) can be done explicitly, with the result that 

&(P )4-i $ _ w = ^[ii|]. 04) 

Here I 3 = \/fSm' and P 2 = Po - Pi- 

We see that both Qn and Q22 give rise to a branch point at p 2 = AM, and become 

complex for p 2 > 4M, whereas they are real for p 2 < AM. Furthermore, Q22 (p) has a 

simple zero in p 2 at p 2 = 0, which means the corresponding Green's function has a pole. 



Together, these results suggest that the fermion acquires a mass mp = vM, while the 
pole at p 2 = is evidence for the existence of a would-be Goldstone boson that reflects 
the condensation of (ijjip)- 

In this note we have analyzed a 2-dimensional model that exhibits the formation of 
Cooper pairs in leading order in i. This condensation occurs for all values of the coupling 
(as long as the bare coupling is positive) for any value of the chemical potential n, including 
H = 0. The coupling itself is asymptotically free. At \i = 0, we have computed the two- 
point functions associated with the composite fields ipip and and have found that tp 
acquires a Majorana mass 2Bq, where So = (V'?/')- We also find evidence for a massless 
pole, which indicates the spontaneous breaking of fermion number at large N. However, 
we expect that the {ipip) condensate will be unstable against higher order corrections in 
jj, so as not to violate the Coleman-Mermin- Wagner theorem. Bearing this in mind, we 
speculate that our GiV-like model at large N could serve as a theoretical laboratory for a 
one-dimensional superconductor [9]. 
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